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ABSTRACT 

While X-ray measurements have so far revealed an increase in the volume-averaged 
baryon fractions fb(r) of galaxy clusters with cluster radii r, fb(r) should asymptot- 
ically reach a universal value /h(oo) = provided that clusters are representative 
of the Universe. In the framework of hydrostatic equilibrium for intracluster gas, we 
have derived the necessary conditions for Jh(oo) = The X-ray surface brightness 
profile described by the (3 model and the temperature profile approximated by the 
polytropic model should satisfy 7 « 2(1 — 1/3(3) and 7 « 1 + 1/3/3 for (3 < 1 and 
/3 > 1, respectively, which sets a stringent limit to the polytropic index: 7 < 4/3. 
In particular, a mildly increasing temperature with radius is required if the observa- 
tionally fitted (3 parameter is in the range 1/3 < (3 < 2/3. It is likely that a reliable 
determination of the universal baryon fraction can be achieved in the small (3 clus- 
ters because the disagreement between the exact and asymptotic baryon fractions for 
clusters with (3 > 2/3 breaks down at rather large radii (<; 30r c ) where hydrostatic 
equilibrium has probably become inapplicable. We further explore how to obtain the 
asymptotic value /b(oo) of baryon fraction from the X-ray measurement made pri- 
marily over the finite central region of a cluster. We demonstrate our method using a 
sample of 19 strong lensing clusters, which enables us to place a useful constraint on 
fb(oo): 0.094 ± 0.035 < /b(oo) < 0.41 ± 0.18, corresponding to a cosmological density 
parameter 0.122 ±0.069 < Sl M < 0.53 ±0.28 for H = 50 km s" 1 Mpc" 1 . An optimal 
estimate of /&(oo) based on three cooling flow clusters with (3 < 1/2 in our lensing 
cluster sample yields (/&(oc)) = 0.142 ± 0.007 or Q M = 0.35 ± 0.09. 

Key words: cosmology: theory — galaxies: clusters: general — intergalactic medium 
— X-rays: galaxies 



1 INTRODUCTION 

In the standard scenario of structure formation, a typical 
galaxy cluster draws its matter (baryon + nonbaryon) from 
a region of radius of ~ 20 Mpc in the Universe. There- 
fore, it is widely believed that clusters should be fair sam- 
ples of baryonic and nonbaryonic matter compositions, and 
thus their baryon fractions ft can be used to determine the 
average mass density of the Universe, Qm, in conjunction 
with the Big Bang Nucleosynthesis (BBN) (e.g. White et 
al. 1993; David, Jones & Forman 1995). However, all the 
X-ray measurements have so far shown an increase in the 
baryon fractions /t,(r) of clusters with radii and no any evi- 
dence for an asymptotic tendency towards a universal value 
at large radii (White & Fabian 1995; Ettori, Fabian & White 
1997; David 1997; White, Jones & Forman 1997; Ettori & 
Fabian 1999; for a recent summary see Wu 1999a). The con- 



flict between X-ray measurement and theoretical expecta- 
tion becomes even more serious when the observed temper- 
ature profiles T(r), which often exhibit a radial decline at 
large radii and can be well approximated by the polytropic 
models of 7 « 1.1-1.3 (Markevitch et al. 1998; Ettori & 
Fabian 1999), instead of an isothermal gas distribution are 
used (Henriksen & Mamon 1994; Henriksen & White 1996; 
Markevitch et al. 1999). In particular, it has been realized 
that the puzzle is unlikely to be associated with the con- 
ventional cluster mass estimates at least within the Abell 
radius, which relies upon the hydrostatic equilibrium hy- 
pothesis for the dynamical state of clusters in the compu- 
tation of their total masses. This point has been recently 
justified by the excellent agreement among the X-ray, op- 
tical and weak lensing determined cluster masses on scales 
0.5 <, r <, 3 Mpc (Allen 1998; Wu et al. 1998), where and 
also hereafter the Hubble constant is taken to be Ho = 50 
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km s _1 Mpc -1 . Yet, we cannot exclude the possibility that 
the regions accessible to current observations are still not 
large enough for the volume-averaged fb(r) in clusters to be 
representative of the Universe. Meanwhile, it remains un- 
clear how accurate the conventional cluster mass estimates 
will be on the outskirts of clusters where the infalling mat- 
ter probably becomes important. A definite resolution to the 
puzzle may require detailed studies of these external regions. 
A good example of such studies has been provided in optical 
by Geller, Diaferio & Kurtz (1999) for the Coma cluster al- 
though their work yields no information about the baryonic 
mass of the cluster. 

In this paper we make no attempt at resolving the 
puzzle. Instead, we first study the necessary conditions for 
/b(oo) = ft if clusters are the well dynamically-relaxed sys- 
tems and share a common value of baryon fraction f b at 
r — > co. We then explore the possibility of deriving the 
asymptotic value /b(oo) from a local measurement fb(r) 
made primarily over the central region of a cluster. Finally, 
we apply our method to a strong lensing cluster sample and 
demonstrate how a useful constraint can be set on the uni- 
versal value fb. 



2 BARYON FRACTION 

2.1 Volume-averaged baryon fraction 

Since the contribution of stellar mass M* to the baryon frac- 
tion of a cluster is typically 5 times smaller than that of gas 
mass M gas (e.g. White et al. 1993; Ettori et al. 1997), we will 
not include M* in the estimate of the cluster baryon fraction 
below. Following the conventional treatment (Cowie, Hen- 
riksen & Mushotzky 1987; Henriksen & Mamon 1994), we as- 
sume a (3 model for intracluster gas characterized by the elec- 
tron number density profile, n e (r) = n e o[l + (r /r c ) 2 ]~ s , and 
an equation of state, T(r) = To[n e (r)/n e o] 7_1 , where n e o 
and To are the corresponding central values and S ^ 3/3/2. 
The total mass in gas within a sphere of radius r is simply 
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ga S (x) = 4nfj, e m p neor c M ga s(x), 
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where x = r/r c , pi e = 2/(l+A), and X — 0.768 is the hydro- 
gen mass fraction in the primordial abundances of hydrogen 
and helium. If the X-ray emitting gas is in hydrostatic equi- 
librium with the underlying gravitational potential of the 
cluster, then the total dynamical mass within r (or x) is 



M tot (r) 
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in which /i; = 0.585 denotes the mean molecular weight. 
The volume-averaged baryon fraction within x is then 



fb(x) = 



M gas (x) 
M tot (x) 



/o(l + ^ 2 ) 1+i(7 - 1) ^#, 



in which 
/o = 



(4) 



(5) 
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We now examine the dependence of baryon fraction on 
cluster radius. In the central region where x — » 0, 



MO) = f • 



(6) 



At large radius where x ^> 1 the baryon fraction varies 
asymptotically as 



f b (x) « /oC^ 1 '" 1 + 



/o 



-X 



2-25+25(7-1) 
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where C is a constant. This disagrees with the previous 
result given by Henriksen & Mamon (1994) who claimed 



fb(x) 



r 2-<5[2-(7-l)] 



at large radius, which results in 7 = 



3 — 2/8 if /b(oo) is asymptotically constant. It is likely that 
they have mistaken x 2+s ^'^ for the asymptotic expansion 
of (1 + x 2 ) 1+s< -~'~ 1 \ and oversimplified the integral M gaa (x). 
For instance, their analysis has neglected the asymptotic be- 
havior of M gas {x) oc In x at large radius for 8 = 3/2. Indeed, 
in the single j3 model with 8 — 3/2 or f3 w 1, /t(co) appears 
to be divergent. For any other 8 (i.e. 8 7^ 3/2), if we require 
/b(oo) to be constant, it is easy to show that the polytropic 
index 7 should satisfy 



7 = 



2-|, 8< 
1 + i S> 



(9) 



(10) 



This sets a stringent limit to the value of 7 
4 

7< 3- 

Namely, the necessary condition for the baryon fraction of a 
cluster asymptotically approaching a universal value at large 
radius is that the intracluster gas has a polytropic index 7 < 
4/3. It appears that all the observed temperature profiles 
of clusters so far have indeed met this simple requirement 
(e.g. Markevitch et al. 1998). Furthermore, because in the 
great majority of cases the observed X-ray surface brightness 
profiles show (3 < 1, we will only focus on the situation of 
where 8 w 3(3/2 < 3/2, for which the relationship between 
the volume-averaged baryon fractions measured at r = 
and r — 00, according to eqs.(6)-(9), is 



Moo) 



3-28 



Mo) 



1 
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Nevertheless, the necessary conditions eq.(9) will yield an 
unphysical value of 7 < if 8 < 1/2, or an unusual result 
of0<7< 1 if 1/2 < 5 < 1. The former may be avoided 
because the observationally fitted (3 parameters from the X- 
ray surface brightness profiles of clusters are usually larger 
than 1/3, i.e., 8 1> 1/2. In the latter case (1/3 < (3 < 2/3), 
the temperature profile of a cluster is required to slightly 
increase with outward radius in order for /b(co) to main- 
tain constant. This last point is indeed beyond our natural 
expectations. Whether or not such a requirement is consis- 
tent with the spectroscopic data will be addressed in the 
discussion section. 

It deserves to examine the issue as to how fast the 
baryon fraction fb(x) approaches the universal value. A rel- 
evant question is: to what radius can we take the volume- 
averaged baryon fraction fb(x) as a good approximation 
of the universal value ? We display in Fig.l the varia- 
tion of fb(x) with radius for 8 = 3/4, 1 and 1.275, re- 
spectively, along with our asymptotic approximations from 
eq.(ll). Surprisingly, the disagreement between the exact 
and asymptotic baryon fractions breaks down at rather 
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larger radii, especially for clusters with large 8. Conse- 
quently, it is unlikely that one can 'directly' measure the cos- 
mic baryon fraction within a cluster of S > 1 (or (3 > 2/3), 
which needs the detection of X-ray emission out to a radius 
of r <; 30r c , corresponding to r <; 7.5 Mpc for an average 
X-ray core radius of 0.25 Mpc. This even does not account 
for the fact that hydrostatic equilibrium fails in this external 
region. Indeed, the problem becomes much more serious for 
larger (3 clusters because the observationally fitted (3 param- 
eter is strongly correlated with r c , with a large value of (3 
giving rise to a large r c value. Moreover, it is apparent from 
Fig.l that fb{x) increases monotonically with radius within 
the region accessible to current observations (r/r c Ss 10), 
which may provide a reasonable explanation for the present 
status of X-ray measurements of cluster baryon fractions 
mentioned at the very onset. We note, however, that it 
should be possible to estimate the universal baryon fraction 
from the small (3 clusters. For example, in a (3 = 1/2 clus- 
ter the deviations of the baryon fractions measured within 
r = 10r c and 5r c from the asymptotic value are only 5.5% 
and 11.5%, respectively. It will be expected that a consis- 
tent baryon fraction can be obtained and applied to the 
determination of the cosmic density parameter Om when an 
ensemble of small (3 clusters are used. 



2.2 Projected baryon fraction 

Next we compute the projected cluster baryon fraction 
fb(b) = m gas (b)/mtot{b) within radius b along the line of 
sight, where m ga s(b) and mtot{b) are the projected gas and 
gravitating mass of the cluster within 6, respectively. Under 
the same working hypothesis as in the above discussion, we 
have 



(x) = 47r/i e m p n e or c m 9as (a;), 
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where x = b/r c . The projected baryon fraction within x is 
then 
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In particular, if we take the condition, 7 = 2—1/5, found 
in the above subsection for S < 3/2 in order for M°°) to 
become a universal value, it can be shown that the above 
expression in the case of 5 > 1/2 (i.e. (3 <; 1/3) reduces to 



Mo) = y- 



(16) 



As a result, in the limits of 1/3 < (3 < 1 and 7 < 4/3, we 
find according to eqs.(6), (11) and (16) 



f , > 2 ~ f (n , „ 2 / t (Q) 
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3 APPLICATION TO STRONG LENSING 
CLUSTERS 

In this section we demonstrate how to constrain the asymp- 
totic baryon fraction of a cluster at large radius using the 
X-ray measurement made within a limiting detection ra- 
dius. For this purpose, we select clusters from the Strong 
Lensing Cluster Sample compiled by Wu et al (1998), which 
contains 48 arclike images of background galaxies gravita- 
tionally lensed by 38 foreground clusters. Here we exclude 
the clusters (l)whose X-ray temperatures are unknown or 
given by indirect methods such as the L x -T and a-T corre- 
lations, where L x is the X-ray luminosity and a is the veloc- 
ity dispersion of cluster galaxies, or (2)whose X-ray surface 
brightness profiles are not well fitted by the (3 models due 
to either the poor data quality (e.g. C12244, etc.) or other 
mechanisms intrinsic to clusters. Nevertheless, we have in- 
cluded the well-known arc-cluster A370 but assumed a core 
radius r c = 0.25 Mpc and (3 = 2/3 although its X-ray sur- 
face brightness is not well constrained (e.g. Ota, Mitsuda 
& Fukazawa 1998; Arnaud & Evrard 1999). This results in 
a sample of 19 clusters, in which there are 25 strongly dis- 
torted images of distant galaxies (Table 1). The projected 
cluster mass interior to the position (r arc ) of an arclike im- 
age can be simply estimated through 



TTllensij' arc) — TTf ' arc^criti 



(18) 



where E crit = (c 2 /4irG)(D s / D d D ds ) is the critical surface 
mass density, with D d , D s and D ds being the angular di- 
ameter distances to the cluster, to the background galaxy, 
and from the cluster to the galaxy, respectively. Spectro- 
scopic data have not been available for about half of the 
arcs, for which we assume a redshift of z s — 0.8. The un- 
certainties of mi en3 due to the unknown redshifts of arcs 
are very minor. While the strong lensing can indeed yield 
the projected gravitating mass of a cluster independently of 
cluster matter contents and their dynamical state, eq.(18) 
may lead to an overestimate of the true cluster mass by a 
factor of ~ 2-4 because it does not take the contribution of 
substructures into account. The factor of ~ 2-4 comes from 
a statistical comparison of cluster masses determined from 
various methods including strong lensing, weak lensing, X- 
ray, optical and numerical simulations (Allen 1998; Wu et al. 
1998; Wu 1999b). Consequently, the ratio of the projected 
mass in gas m sas (r arc ) to the strong lensing derived mass 
mi ens (r arc ) provides a low limit to the true baryon fraction 
fb(r a rc) within r arc . Moreover, since all the arclike images 
are detected within the central regions of the clusters, we 
can use /t(r arc ) to approximately represent the projected 
baryon fractions at cluster centers, /i,(0) £s /b(r arc ). That 
is, 



M0) « Mr a 



s (j*arc) 



(19) 
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Table 1. Strong Lensing Cluster Sample 



cluster Zduster T (kcV) L\ r c c /3 z arc r c arc m d gaa ™f ens m gas /m len3 (%) / (%) e 



A370 





.373 


7 i o+1.05 
' - id -1.05 


20.8 


+2.0 

-2.0 





.25 





.67 


1. 


.3 





.35 


2, 


.811° 


25 
23 


130 


2 16+ - 19 
Z - 1D -0.18 


18.6±|-_? 






















0. 


724 


0. 


16 


0, 


.731° 


.07 
.06 


29.0 


9 Kr, + 0.23 

Z - OZ -0.21 




A963 





.206 


6 1S+ 045 


16.1 


+ 1.8 

-1.8 





.08 





.67 






0. 


0517 


0, 


12l° 


.01 
.01 


2.5 


4 91+0- 33 
^■ y± -0.36 


11 1+ 1 ' 4 






















0. 


711 


0. 


080 


0, 


.271° 


(12 

.02 


6.0 


4 44+0-30 
*-**-0.33 




A1689 





181 




50.5 


+ 1.1 

-1.1 





.13 





.65 






0. 


183 


1, 


.561° 


.03 
.03 


36 


4 09+O.O8 
^•^-O.OO 


15.01°;* 


A1835 





252 


9.80li;l° 


105±\i 





.074 





65 









150 


1, 


651° 


18 
.17 


19.8 


o oc+0.89 
°- oo -0.83 


14.81 4 ;? 


A2163 





.203 


1 4 fi0+°' 85 
rt.DU_ 85 


133±iS 





.34 





.64 


0. 


728 


0. 


0661 





.241° 


.02 
02 


4.1 


5 84+0-50 
°-°^-0.51 


21.4l2 g 


A2218 





.171 


7 10 + ' 20 

' - 1 -0.20 


18.1 


+0.75 
-0.75 





.23 





65 


1. 


034 


0, 


26 


1, 


■61±g 


.04 
.04 


27 


5 98+° 17 
°- ao -0.17 


16 l +0 ' 9 
lu - -0.9 

























702 


0, 


0794 


0, 


181° 


01 

.01 


6.23 


9 qo+0.08 
a -0.08 
























0. 


.515 


0, 


0848 


0, 


■21±8 


01 

.01 


5.70 


3.631°;!° 




A2219 





.228 


rz.4u_ 5() 


64.6 


+7.0 
-7.0 





.15 





,64 









079 





.351° 


.02 
02 


5.6 


6-82l°;« 


12 1+ 1 ' 3 
iZ -l-1.2 


























0. 


110 


0, 


.661° 


1)4 

04 


16.0 


4 41+0-26 
^'^ -0.27 




A OQOfl 

AzoyU 


U. 


.ZZo 


11 - 10 -1.00 


63.5 


+ 14.9 
-14.9 


n 

u. 


no 


u 


aft 

DO 


n 

u. 


Ol Q 

y lo 


U. 


1 77 
111 


1 

1. 


■«„ 


20 
.21 


OK A 
Z0.4 


c: K4 +0.78 
5 - b4 -0.81 


17 ,+4.3 


A2744 





.308 


11 00+ ' 50 

il - uu -0.50 


64.4 


+ 14.4 
-14.4 





.45 





.65 






0. 


1196 


0, 


•sit? 


Ofi 
.07 


11.36 


4 4»+0- 53 


24.81 4 ;? 


AC 114 





310 


9 76+ 1 ' 04 

a -'°-0.85 


38.1 







.30 





67 


1, 


86 


0, 


0676 


3 


47+0 

-0 


.08 
.0!) 


2.98 


9 C7+O.O6 
z - D '-0.07 


18.7ll;f 

























.639 





.35 





16+° 

±u -0 


.00 
00 


130 


5 40+0-12 
°' 4Z -0.14 




C10500 





.327 


7 20+ 3 ' 70 
'• ZU -1.80 


17.5 


+ 1.3 
-0.9 





.030 





.36 






0. 


.15 


1, 


■21±S 


14 
.15 


19.0 


6 S7+ ' 07 


1 0+0.6 

l- J -0.5 


MS0440 





.197 


5 30+ 1 ' 27 
°- JU -0.85 


7 4+1-0 

'•^-1.0 





.054 





68 


0, 


530 





.089 





231° 


03 
.03 


8.9 


9 Kfl+0 29 
2 - 60 -0.31 


j r+2.5 
'•°-2.2 


MS0451 





.539 


10.1711;" 


53.7 







.088 


0. 


.84 






0. 


190 


1, 


sit? 


.05 
.05 


52 


2 91+0-1° 
z ' a -0.10 


14.011J 


MS1008 





306 


7.29+?'« 

— 1.52 


9.1± 


1.2 
1.2 





12 


0. 


.68 






0, 


26 


1, 


18tS 


15 
.16 


61 


1.94±°;» 




MS1358 





.329 


6 5 +0 - 7 


21.8 


+3.8 
-3.8 





.023 





36 


4 


92 


0. 


121 


1. 


.021° 


11 
.12 


8.27 


12.411-3 


l-4t° j 


MS1455 





257 


c 4 r; + 0.29 


29.4 


+ 1.5 
-1.5 





.040 





66 






0, 


098 


0. 


521° 


(12 
.02 


8.6 


6 04+ ' 24 
°- u4 -0.23 


I2.0ll;g 


MS2137 





313 


4.371°,;?* 


26.3 


+3.6 
-3.6 





.054 





,64 









0874 


0. 


44+0 


05 

04 


7.1 


6 17+°- 7 i 

O,1 '-0.56 


16.6+5° 


PKS0745 





.103 


8.7±l:| 


57.2 


+4.4 
-4.4 





.056 





,59 


0, 


433 


0. 


0459 


0. 


171° 


01 
01 


3.0 


5-8itS:S 


9-4111 


RXJ1347 





.451 


11.37l ; 9 ° 


1981H 





.024 





.42 









24 


6, 


.151° 


47 
.48 


42 


14.611- 


3.4l°; 6 6 



"Lensing and X-ray data (except for r c and 0) are taken from Wu et al. (1998) and Wu, Xue & Fang (1999), respectively; 
b (X-ray) bolometric luminosity in units of 10 44 erg s — 1 ; 
c in units of Mpc; 

^projected mass within r arc in units of 1O 13 M0; 
e calculated by eq.(5). 



The resultant values of m gas (r arc ), m( ens (r arc ) and their 
ratios within the positions of 25 arcs have been given in 
Table 1 for a flat cosmological model Q. — 1. 

On the other hand, we can also calculate the volume- 
averaged baryon fraction around r — for each lensing 
cluster using eqs.(5) and (6), along with the restriction 
7 = 2 — 1/5. Unfortunately, the temperature profiles have 
not been obtained for these lensing clusters. Instead, present 
X-ray observations have only provided the global emission 
weighted-temperatures 



J a(T)T(r)n 2 e (r)r 2 dr 
f a{T)nl{r)r 2 dr ' 



(20) 



where a(T) is the cooling function. When T is used to calcu- 
late /o, we will either overestimate (for f3 > 2/3) or underes- 
timate (for P < 2/3) the values of fo because T are smaller 
(greater) than the central values To in polytropic models 
with (3 > 2/3 (/3 < 2/3). We have attempted to assign the 
upper and low limits to the central baryon fractions, respec- 
tively, for the 6 and 13 clusters with /3 > 2/3 and f3 < 2/3 
in Table 1 according to 



A(o) = -/o| To= 



(21) 



However, it appears that we cannot make a clear distinc- 
tion between the resultant upper and low limits among the 
19 clusters (see also Fig. 2). This probably implies that the 
evaluation of fb(0) is little affected by the approximation of 
the emission weighted temperature to the central tempera- 
ture. 

The low limits fb(0) [eq.(19)] and the central baryon 
fractions fb(0) approximated by eq.(21) can be converted 
into the total baryon fractions of clusters /&(co) in terms of 
the relations established in the above section [eqs.(ll) and 
(17)], namely, 



fb,lens < fb (OO) « /(,,; 



where 

fb,lens 



TFlgas arc) 



fb.xray — 



3(1-/3) TFllens {Tare) 

1 



3(1-/3) 



/ol 



T =T- 



(22) 

(23) 
(24) 



In Fig. 2 we plot the results of fb,iens and fb,xra V obtained 
from the 25 arclike images among 19 lensing clusters listed 
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in Table 1. Note that fb,xray is independent of arc position 
r a rc, and here we utilize r arc to represent the lensing cluster 
only. It is apparent that the lensing and X-ray results are 
clearly separated except for three clusters, C10500, MS1358 
and RXJ1347. Dispersion in fb,iens is relatively small, and 
averaging over the 25 data points yields 



Moo) > </ Mens ) = 0.094 ± 0.035, 



(25) 



where the error bar has not accounted for the uncertainties 
of (3 and r c arising from the (3 model fitting. The average 
value (fb,iens) increases only slightly if the cooling flow clus- 
ters alone are considered: 



{Mens) =0.106 ± 0.033. 



(26) 



The reason why we emphasize this result is that an excellent 
agreement between cluster masses determined from strong 
lensing and X-ray measurements for cooling flow clusters has 
been claimed by Allen (1998). If this is true, we would expect 
Mens = M°°)- Nevertheless, we will still adopt eq.(25) 
for caution's sake: First, the recent study of Lewis et al. 
(1999) on the cluster mass estimates using the CNOC cluster 
sample does not reveal the large discrepancy between cooling 
and non-cooling flow clusters; Second, the values of fb,iens 
we have found for cooling and non-cooling flow clusters are 
essentially consistent with each other (see Fig.2). 

A glimpse of Fig.2 reveals that the three clusters 
(C10500, MS1358 and RXJ1347) which give smaller values 
of fb.xray all have their (3 values less than 1/2. This mo- 
tivates us to examine the dependence of fb,x ray upon /3 
for our cluster sample (Fig. 3). Indeed, the values of fb,x ray 
derived from the volume-averaged baryon fractions among 
clusters with (3 > 2/3 are systematically high. As has been 
pointed out in section 2.1, these asymptotic values cannot be 
used as a reliable indicator of the universal baryon fraction 
because the agreement between the exact and asymptotic 
values of the baryon fractions occurs at rather large radii 
r <; 30r c where the hydrostatic equilibrium hypothesis may 
have already broken down. On the other hand, there are 



good reasons that the asymptotic values f b ^ 



obtained 



among C10500, MS1358 and RXJ1347 can be considered as 
a good approximation of the universal baryon fraction: First, 
the (3 parameters for these three clusters are smaller than 
1/2, which guarantees the condition M x ) ~ fb,xray within 
~ 10r c ; Second, all these three clusters are classified as the 
cooling flow clusters, for which the total cluster masses M to t 
could be determined relatively accurately from eq.(3) (Allen 
1998); Third, the asymptotic values fb,x ray from the volume- 
averaged baryon fractions among these three clusters are es- 
sentially consistent with the asymptotic values fb,iens , and 



(Me™) =0.121 ±0.088; 

(Mray) = 0. 142 ± 0.007. 



(27) 
(28) 



As a conservative estimate, we utilize the average values 
(fb,iens) and (fb,xra y ) among all the clusters to be the low 
and upper limits to the universal baryon fraction 



0.094 ± 0.035 < / 6 (oo) < 0.41 ± 0.18, 



(29) 



which, in combination with the BBN prediction (Walker et 
al. 1991), corresponds to the following constraints on the 
cosmological density parameter 



where the uncertainty of the BBN prediction, AQb = 0.01, 
has been included, and Qb is the average baryon mass density 
of the Universe in units of the critical density for closure. If 
the average value (fb,xr ay ) in the three cooling clusters with 
(3 < 1/2 is adopted, we have 



tl M =0.35 ±0.09. 



(31) 



0.122 ± 0.069 < Q.m < 0.53 ± 0.28, 



(30) 



Finally, we investigate whether the asymptotic baryon 
fractions (or limits) fb,iens and fb,x ray derived from the 
strong lensing clusters depend on X-ray temperature. This 
is important because variation of baryon fractions among 
clusters of different temperatures is not allowed in the stan- 
dard models of structure formation. For the time being we 
can only take a less serious approach to the problem in the 
sense that the exact baryon fractions of these lensing clusters 
have remained unknown. Fig.4 illustrates the dependence of 
fb,iens and fb,xra y upon the emission-weighted temperature 
for the 19 clusters in Table 1. The best-fitted fb,i en s — T and 
fb, X ray — T relations are, respectively, 

f 1n -1.38±0.267p0-39±0.28 , . 

Jb,lens = 1U 1 ; (62) 

p 1r ,-0.87±0. 357^0. 47±0. 38 . . 

Jb,xray = 1U 1 , (66) 

where T is in units of keV. These relations are marginally 
consistent with no increase of baryon fractions with T as 
predicted in the standard models although they cannot be 
used as a direct evidence. Similar result has recently been 
reported by Mohr, Mathiesen & Evrard (1999). 



4 DISCUSSION AND CONCLUSIONS 

The baryon fraction fb{r) of a cluster should be averaged 
over a sufficiently large cluster volume in order for fb(r) to 
be representative of the universal value /&. However, this 
is limited by the sensitivity of current detectors which can 
only probe a finite region of the X-ray cluster. We have thus 
studied the condition for fb = fbipo) and also the possibility 
of deriving /b(oo) from the central baryon fraction fb(0) or 
fb{0)- It has been shown that the polytropic index of intra- 
cluster gas should satisfy 7 < | if the baryon fraction of a 
cluster is required to asymptotically reach a universal value 
fb at large radius. It appears that all the clusters whose tem- 
perature profiles have been available thus far indeed meet 
this simple requirement. We have found that clusters with 
small (3 parameters are likely to provide a better estimate 
of the universal baryon fraction, which arises because the 
exact and asymptotic forms of the baryon fractions merge 
at rather large radii r i> 30r c for f3 <; 2/3, where hydro- 
static equilibrium becomes to be questionable. We have also 
demonstrated how to set useful constraints on the baryon 
fraction fb , and thereby the average mass density of the Uni- 
verse, fljf, using a sample of 19 strong lensing clusters. Our 
analysis gives support to a medium density universe, and 
the scatters in the present estimate of Q.m can be greatly 
reduced when a large cluster sample is employed. 

However, there are several factors which should be taken 
seriously. Our major concern is the current working model 
for clusters, in which the hydrostatic equilibrium hypoth- 
esis and the conventional f3 model for intracluster gas are 
extrapolated to sufficiently large radii. The availability of 
this working model to the Abell radius has been demon- 
strated by the reasonable agreement between different mass 



© 0000 RAS, MNRAS 000, 000-000 



6 Wu and Xue 



estimators. In principle, clusters can approximately be re- 
garded as dynamically-relaxed systems out to the falling 
shock radii where cluster matter is mixed with the infalling 
matter from background universe. Presumably, our conclu- 
sion cannot be applied beyond these radii, within which the 
volume-averaged matter composition should be of cosmo- 
logical significance. Even so, since these boundaries are well 
beyond the regions accessible to present observations, it is 
unclear whether the intracluster gas in the outmost regions 
still follows the j3 and polytropic models. 

To ensure the universal constancy of the cluster baryon 
fraction at large radius, the intracluster gas should have a 
polytropic index of < 7 < 1 if the fitted (3 parameter to 
the X-ray surface brightness of the cluster is in the range 
of from 1/3 to 2/3. This requires an increasing tempera- 
ture with cluster radius, in conflict with the naive specula- 
tion that the overall temperature profile of a cluster should 
drop with radius. The latter occurs, according to our pre- 
diction, only inside cluster with 2/3 < j3 < 1. At present, 
the spectroscopic analyses have resulted in a discrepancy in 
the radial temperature gradients due primarily to the prob- 
lem with the sensitivities and resolutions of detectors. For 
instance, by contrast to the remarkably radial temperature 
decline claimed by Markevitch et al. (1998), the nearly flat 
temperature profiles have recently been reported by Irwin, 
Bregman & Evrard (1999) based on a detailed comparison 
of the ASCA and ROSAT PSPC determined temperature 
profiles of 26 clusters. In particular, the trend of an increas- 
ing temperature with radius is likely to present in a few 
cases (e.g. 2A0335, A401, A1651, A3571, etc.). Therefore, it 
is not impossible that the temperature profiles would have a 
mild increase with radii for clusters with 1/3 < f3 < 2/3. A 
conclusive test for this claim will be provided by the space 
experiments like AXAF and XMM. 

Another concern is the poorly determined (3 and r c pa- 
rameters used in the computations of the baryon fraction 
limits fb,iens and fb,xray, because /(,(oo) depends sensitively 
on these two parameters. To date, the X-ray surface bright- 
ness profiles of clusters have not been well measured at large 
radii (<; 1 Mpc). Therefore, there may exist rather large 
uncertainties in the present fitting of /3 models, apart from 
the fact that different observations at different energy bands 
might give very different values of (5 and r c even for the same 
cluster. In particular, the f3 model fitting is affected by the 
procedure of whether or not the cooling flow regions are rea- 
sonably excluded. For instance, from a detailed study of the 
surface brightness distributions of 25 cooling flow clusters, 
Vikhlinin, Forman & Jones (1999) have shown that the j3 
values are systematically higher outside the cooling flow re- 
gions. Alternatively, a double (5 model has been often used 
in recent years to represent quantitatively the excess emis- 
sion in the central cores of cooling flow clusters (e.g. Ikebe et 
al. 1996; Xu et al. 1998; Mohr, Mathiesen & Evrard 1999), 
which also results in a larger value of /3 for the extended gas 
component. So, it is possible to release the requirement that 
X-ray temperature should increase with radius in order for 
fb(oo) to remain constant if large j3 parameter is adopted. 
Furthermore, because our attempt to derive the universal 
value /f,(oo) from the central baryon fractions /t.(0) depends 
sensitively on the central gas distributions, our current re- 
sults are subject to the presence of the sharp peaks in emis- 
sion concentrated in the cores of some clusters. 



Nevertheless, if clusters are representative of the Uni- 
verse, and if the entire intracluster gas is in hydrostatic equi- 
librium with the underlying gravitational potential of clus- 
ters and follows the (5 model, we would expect that the X- 
ray surface brightness and temperature profiles of clusters 
should satisfy the necessary condition eq.(9), which reads 
7 w 2(1 - 1/3/3) for 1/3 < fi < 1. This will be directly 
testable when these two profiles are well measured (e.g. by 
AXAF, XMM, etc.). Even if the temperature profile of a 
cluster is not available, we will still be able to obtain its 
baryon fraction according to eq.(ll) or (17) by accurately 
measuring its X-ray surface brightness distribution and its 
central temperature To, especially for a cluster with smaller 
(3 parameter. Applications to existing data and future ob- 
servations will be made in our subsequent work. 
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Figure 2. Asymptotic baryon fractions of 19 lensing clusters in Table 1 plotted against arc positions r arc . Filled symbols are the low 
limits derived from the projected baryon fractions within the 21 arc positions, and open symbols are the X-ray results derived from /b(0). 
For the latter, r arc are used to denote the clusters only (we take the first arc position for the multiple-arc system). We use the squares 
to denote cooling clusters and the triangles, non-cooling flow ones. Dotted and dashed lines represent the average values of the lensing 
and X-ray results, respectively: {fb,iens) = 0.094 ± 0.035 and {fb,xray) = 0.41 ± 0.18. 
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Figure 3. Dependence of asymptotic baryon fractions fb, X ray upon (3 parameters for 19 lensing clusters. 
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Figure 4. Asymptotic baryon fractions of clusters vs. X-ray temperature for 19 lensing clusters. Filled circles are the low limits from 
gravitational lensing, and open squares are the X-ray results. Dotted and dashed lines are the best-fitted fb~T relations to the two sets 
of data, respectively. 



© 0000 RAS, MNRAS 000, 000-000 



